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ABSTRACT 

The flow pattern and thermal stratification of a cryogenic cylindrical 
tank are numerically studied. The tank sidewall is subjected to either a uni- 
form heat-flux or two discrete levels of uniform heat-flux at the upper and 
lower halves of the tank wall. The tank bottom is kept at a constant tempera- 
ture controlled by the heat exchanger of a thermodynamic vent system. The tank 
pressure is also assumed constant resulting in a constant saturation tempera- 
ture at the interface which is higher than the tank bottom temperature. The 
effects of vapor motion and vapor superheat on the mass and heat transfer proc- 
esses at the interface are assumed negligible such that the calculations of 
liquid region can be decoupled from the vapor region. Dimensionless steady- 
state conservation equations are solved by a finite-difference method. The 
effects of modified Rayleigh number, Prandtl number, tank aspect ratio, wall 
heat-flux parameter, and wall heat-flux distribution on the liquid velocity and 
temperature fields are investigated. Also, their effects on the rate of heat 
transfer through the interface and the tank bottom are examined. 


NOMENCLATURE 

A surface area of the interface (cross-section area of the tank) 

Ar aspect ratio, H/D 

Bo tank Bond number, pgD 2 /4o 

Cp specific heat at constant pressure 

D tank diameter 

Gr* modified Grashof number, p 2 (igq w H 4 /kp 2 
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Pr 
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Ra* 

Rq 
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gravitational acceleration 

gravitational acceleration at the Earth surface 

liquid level height 

latent heat of evaporation 

Jakob number, C p (T s - Tt,)/hf g 

thermal conductivity 

evaporation mass flux 

wall heat-flux parameter, q w D/k(T s - T^) 

defined as q w iD/k(T s - Tj,) 

defined as q w u D/k(T s - Tjj) 

local Nusselt numbers at the interface and tank bottom, 
respectively 

pressure 

equilibrium hydrostatic pressure 
dimensionless pressure, (p - pg)/p(a/D ) 2 
Prandtl number, Cpp/k 

heat fluxes out of the liquid through the interface and tank 
bottom, respectively 

wall heat-flux for uniform heating or average wall heat-flux 

wall heat fluxes at the upper and lower halves of the tank wall, 
respectively 

radial coordinate measured from the centerline 
dimensionless radial coordinate, r/D 
modified Rayleigh number, Gr*Pr 

heat flux distribution parameter (the ratio of wall heat-flux 
distribution), q WiU /q W) i 

temperature 

interface temperature 

tank bottom temperature 
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T* 
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u* 

V 

v s 

V* 

V* 
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X* 
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dimensionless temperature, (T - ) / (T s - Tb) 

axial velocity 

dimensionless axial velocity, u/(a/D) 
radial velocity 

radial velocity at the interface 

dimensionless radial velocity, v/(a/D) 

dimensionless radial velocity at interface 

axial coordinate measured from the tank bottom 

dimensionless axial coordinate, x/D 

thermal diffusivity, k/pC p 

coefficient of thermal expansion 

dynamic viscosity 

liquid density 

surface tension 

dimensionless stream function 

average value over the interface or tank bottom 


Subscripts 

b evaluated at tank bottom 

1 evaluated at the lower half of the tank wall 

s or t evaluated at interface (top surface) 

u evaluated at the upper half of the tank wall 


INTRODUCTION 


Space missions in the coming decades will involve the storage and transfer 
of large volumes of cryogenic liquids such as hydrogen under orbital condi- 
tions. The current interest in cryogenic fluids storage and transfer stems 
part from NASA's plan for an ambitious SEI (Space Exploration Initiative) 
including manned voyages to the moon and Mars. One of the key problems in 
storing cryogenic liquids under low-gravity conditions is the pressure control 
of storage tanks without significant loss of the fluid. 


The thermal environment 
insulated cryogenic storage 


in space can result in heat transfer through an 
tank causing fluid temperature stratification and 
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interfacial heat and mass transfer. As a consequence, the tank pressure is 
continually increasing. The prediction of the tank pressure rise rate and the 
technology development for the tank pressure and temperature control systems 
require knowledge of the temperature distribution in the liquid region. Gen- 
erally, cryogenic storage systems of practical interest are dimensionally 
large. Therefore, the effect of natural convection on the thermal stratifica- 
tion process may be significant even in a micro-gravity environment. As for 
example, for a liquid hydrogen storage system, the modified Rayleigh number may 
vary in the range of 10 to 10^ in an accelerational environment of 10 - ® g e or 
less. The problem described in this paper simulates some aspects of a cryo- 
genic fluid storage system with a pressure control device such as a thermody- 
namic vent system (ref. 1). 

The use of a thermodynamic vent system (TVS) is one of the major technolo- 
gies (ref. 2) being developed by NASA for controlling cryogenic tank pressure. 
In a thermodynamic vent system, a small amount of cryogenic liquid is continu- 
ously withdrawn from the tank and passed through a Joule-Thomson valve result- 
ing in a lower pressure and temperature two-phase fluid. The mixture is then 
introduced into a heat exchanger to either cool the tank liquid by absorbing 
its thermal energy and/or to intercept the heat leak which will otherwise be 
transferred into the liquid through the tank walls. The two-phase mixture is 
continually evaporating and the resulting vapor is then vented overboard. 
Therefore, if the TVS heat exchanger is mounted on some portion of the tank 
wall, the temperature of that portion of the tank wall will be relatively cold 
and constant because of the two-phase mixture of the TVS flow. 

Generally, if the temperature of the lower portion of the tank is higher 
than that of the upper portion of the tank, buoyancy- induced mixing will be 
enhanced and the maximum temperature of the fluid is decreased. In this study, 
the tank pressure is assumed constant corresponding to a constant saturation 
temperature at the liquid-vapor interface. The TVS heat exchanger is assumed 
to be mounted on the tank bottom such that the tank bottom temperature is lower 
than the interface temperature. This may represent the worst situation in 
which the maximum fluid temperature is increased. The tank sidewall is subjec- 
ted to a prescribed heat flux. The heated liquid near the tank sidewall flows 
toward the interface under the action of the buoyancy force. A layer of highly 
stratified liquid is then formed underneath the interface at a temperature 
higher than the bulk liquid. At steady state, the energy coming through the 
tank sidewall must equal the total energy transferred out through the interface 
and the tank bottom. The objective of the present study is to numerically 
investigate the effects of buoyancy force, tank aspect ratio, wall heat-flux, 
and fluid thermodynamic properties on the rate of heat transfer through the 
free surface and the cold tank bottom. Their effects on the liquid velocity 
and the temperature fields are also examined. 


PROBLEM FORMULATION 

The physical system considered and coordinates are illustrated in 
figure 1. A circular cylindrical tank of diameter D is filled with a cryo- 
gen such as liquid hydrogen, nitrogen, or oxygen at a liquid level height H. 
The tank sidewall is subjected to either a uniform heat-flux (q w ) or two dis- 
crete levels of uniform heat-flux (q WjU ) and (q Wi i) at the upper and lower 
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halves, respectively, of the tank wall such that the average .all heat-flux 
( n + a i)/2 is equal to q*. A parameter Rq defined as qw.u^w.l 
be W uSed ^characterize the effect of wall heat-irfux distribution The tank 
bottom is maintained at a cold constant temperature Tb controlled by the TVS 
heat exchanger. It is assumed that the effect of vapor superheat and vapor 
motion on the liquid-vapor interface heat and mass transfer is negligib e. 

Thus, the calculations of the liquid region can be decoupled from the vapor 
region. The interface is assumed to be flat, wave free, and at a constant 
temperature T s corresponding to the constant tank pressure. It is fur er 
assumed that the change of liquid level due to the liquid consumption of the 
TVS and mass transfer at the interface is negligible. 

Under the conditions and assumptions described above, the problem essen- 
tially reduces to a steady-state natural convection flow in an enclosure with 
zero friction at the free surface. The dimensionless form of constant-property 
continuity, momentum, and energy equations governing this axisymmetric laminar 
natural convection flow with Boussinesq approximation can be written as 
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Note that the gravity term has been substracted from the momentum equations by 
using the hydrostatic equilibrium equations. The parameters appearing in the 
governing equations are the Prandtl number (Pr), the modified Rayleigh number 
(Ra*), the wall heat-flux parameter (Nh) , and the tank aspect ratio (Ar). 


C u 

Pr = ^- = -j~, 
pa k 


Ra* = Gr*Pr = 


p 2 |3gq w H 4 C p 
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Boundary conditions are required to solve the governing partial differential 
equations (1) to (4). At the centerline, the symmetric conditions are applied: 


Nonslip conditions are used for the tank sidewall and bottom: 


u* = v* = 0 (8) 

The thermal condition at the tank sidewall is either, for a uniform wall heat 
flux 


3T* 

gpr = Nh at the entire tank sidewall 


( 9 ) 


or, for two discrete levels of uniform wall heat fluxes, 


9T* 

3r* 


Nh u at the upper half of the tank sidewall 
Nhj at the lower half of the tank sidewall 


( 10 ) 


where 


Nh u + Nhj = 2 Nh 


(U) 


The liquid vapor interface and tank bottom are maintained at constant tempera 
tures T s and T^ , respectively, with T s > T^. Thus, 


T* = 1 at 



( 12 ) 


and 


T* = 0 at x* = 0 ( 13 ) 

The evaporation- induced (or may be condensation- induced) axial velocity at the 
interface is calculated by 


u 


e 




(14) 


where m e is the interface evaporation mass flux and hf„ is the latent heat 
of evaporation. The nondimensional evaporat ion- induced velocity u* can be 
expressed as e 


u 


* 


e 



(15) 
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where Ja is the Jakob number defined by C p (T s - T b )/hf g . In practica 
applications, the value of Ja for liquid hydrogen, nitrogen, or oxygen, is 
small (<0 . 1 ) and u* is generally negligible. Thus, it is assumed u ? 
fi.e., Ja = 0) in the present calculations. The shear-free condition is a so 
applied at the interface since the vapor viscosity is usually much smaller 
the liquid viscosity. Therefore, the following boundary conditions, 


than 


u* = 0, 


9v* 

3x* 


0 at 


* _ H 

x - D 


(16) 


are used at the interface. 

The energy balance in the liquid region can be expressed as 


J "t “s + I 


% dA b 


q irDH = 0 


(17) 


where q w and qb are the incoming heat flux through the tank sidewall and 
the heat flux transferred out through the tank bottom, respectively, and <jt 
is the heat flux transferred through the interface. The va ues of % 
q h are always taken positive. However, the value of qt t is .^ ak ® n n ^ 
and negative for heat transferred out of the liquid and into the liquid, 
respectively. Equation (17) simply shows that the nej_heat a \ the 
liquid region is zero. The average Nusselt numbers, Nu t and Nu b , at the 
tank interface and tank bottom, respectively, are defined as 


Nu 


t 



(18) 


and 



(19) 


where qt and qh are the average values of q t ,and q b over surfaces 
of the iJterface Snd tank bottom, respectively. In equation 18), the posit ve 
sign is used if q t > 0 and the negative sign is used if q t < 0. By using the 
definitions of Nh and equations (18) and (19), equation (17) 


then becomes 


4 Nh 

Ar - Nu t 

if 

^t 

> 0 

Nu, = { 
b ] 4 Nh 

Ar + Nu 

if 

q t 

< 0 

t 


t 



( 20 ) 


The ratio, Nu b /(4 Nh Ar), represents the fraction of the incoming energy trans 
f erred througLthe cold’surface (i.e., removed by the TVS heat exchanger). If 
the value of Nub/(4 Nh Ar) is greater, then the energy transferred to the 
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interface to cause liquid vaporization is less. The energy ratio can be 
expressed as a function of the relevant parameters 


Nui 


4 Nh Ar = f(Ra *’ Pr - Nh ’ Ar - R q ) 


( 21 ) 


The effects of the above parameters on the ratio, Nu,/4 Nh Ar, and the velocity 
and temperature fields are numerically investigated. 


COMPUTATIONAL METHOD 

f overnin 8 equations (1) to (4) are solved by a finite-difference 
methocL A staggered grid system is employed with the scalar properties, p* 

* ’.*° catet * mi dway between the u* and v* velocity grid nodes. The 
t ini te-di f ference equations are then derived by integrating the differential 
equations over an elementary control volume surrounding a grid node appropriate 
ror each dependent variable. Bounded skew hybrid differencing (BSHD) (ref. 3) 
is used in the convective terms and the integrated source terms are linearized 
Both of these practices enhance numerical stability. To accelerate conver- 
f re obtained by using the pressure implicit split operator 
(P1S0) method (ref. 4) to solve a pressure-correction equation which yields the 
pressure change needed to acquire velocity changes to satisfy mass continuity 
the finite-difference equations are iteratively solved by the alternating 
direction implicit (ADI) method (ref. 5) with under-relaxation until the solu- 
tions are converged. 

A nonuniform grid distribution is generated by using Roberts' transforma- 
tion (ref. 6) with concentration of the grid nodes in the near-wall, near- 
bottom, and near-interface regions where the gradients of flow properties are 
large Several grid distributions (36 by 20, 48 by 27, 60 by 34, and 72 by 41) 
*u 6n ^ rie( ^ t0 °btain reasonable grid- independent solutions. It is found 
tliat the changes in the average radial velocity at the interface (v*) and the 

average Nusse It number at the tank bottom (Nu b ) between grid distributions 

?? 3n< ! 7 \ by 41 are ‘O' 1 P ercen t for the worst case of uniformly-heated 

ank (Nh 4, Ar - 1.25, and Ra* = 10 9 ). Thus, it is decided to use 60 by 34 
grid distribution for all of the calculations in the present study. 


RESULTS AND DISCUSSION 


Numerical solutions are obtained for the aspect ratio (Ar), the modified 
Rayleigh number (Ra ), and the wall heat-flux parameter (Nh) ranging from 
0.25 to 1.25, 10S to l 0 9 t and 1 to 4, respectively. The heat-flux distribution 
parameter (Rq) is varied from 0 to 3. The Prandtl number (Pr) is varied from 

1.25 to 2.6 fo cover the typical values of Pr for liquid hydrogen, nitrogen, 
and oxygen. 6 


Effects of Ra* and Pr 

The effects of the modified Rayleigh number (Ra*) on the flow pattern are 
shown in figures 2(a) to (c) for a uniformly-heated tank with Ar = 1 , Nh = 4, 

8 



and Pr = 1.25. The dimensionless stream function, y, in the figure is calcu- 
lated from 


i|/(x\ r* = 0) = x|/(x* = 0, r* 


X* 

0) - [ r*v* dx* 
J 0 


Ui(x* , r*) = y(x\ r* = 0) + [ r*u* dr* 

J 0 


( 22 ) 


(23) 


with H*( x * = 0, r* = 0) set to zero. For Ra* = 10^ (fig. 2(a)), the small 
velocity field forms a big and relatively symmetric counter-clockwise vortex 
flow pattern and the heat transfer is dominated by both convection and con- 
duction. As Ra* increases, a secondary flow of small clockwise vortex is 
generated in the lower left corner (fig. 2(b)) forcing the center of the 
counter-clockwise vortex moving toward the upper right corner and resulting in 
a narrower area for the boundary layer flow to turn radially toward the centra 
region. The fluid velocity in the boundary and near the interface region is 
greater for larger Ra*. Even though the radial flow turning area near the 
interface is smaller for higher Ra* , the returning flow velocity at the cen- 
terline is still larger than that for lower Ra*. The clockwise vortex is 
lifted detaching from the tank bottom as Ra* further increases and a thin 
boundary layer along the tank sidewall is developed (fig. 2(c)). 

Without buoyancy force (Ra* = 0), there will be no fluid motion and the 
temperature distribution will result from the heat conduction only. As Ra 
increases, the natural covection process due to buoyancy results in warm flui 
moving upward along the tank sidewall which increases the rate of heat trans- 
fer out of the interface. However, an increase in Ra* convects the heat 
from the near-wall region faster. This results in lower and more uniform liq~ 
uid temperature in the near-interface region (as seen from figs. 2 and 3) which 
decreases heat transfer rate_at the interface. These two compensating ef sets 
result in the energy ratio, Nu^/(4 Ar Nh) , only slightly changing with Ra 
for Ra* = 10 5 to 10 9 (as seen in fig. 4). 

Figures 3(a) to (c) show isotherm plots (T*) for the cases corresponding 
to figure 2. Since buoyancy- induced flow for higher Ra* convects heat from 
the tank wall faster than lower Ra*, the maximum temperature is lower for 
higher Ra* and is located at the tank sidewall closer to the free surface. 

A natural convection boundary layer is clearly formed for higher Ra w it 
slight temperature variation outside of the boundary layer. Most of the region 
has temperature higher than the interface. The higher Ra* yields relatively 
more uniform liquid temperature and lower degree_of superheat in the liquid 
region. Figure 4 shows the effect of Ra* on Nuj)/(4 Ar Nh) , which repre- 
sents the fraction of the total incoming energy transferred through the cold 
surface (tank bottom), for a uniform wall heat-flux of Nh = 4, Pr = 1 - 25 < 
and Ar = 1 . The energy ratio reduces from 0.44 for Ra* = 10° to 0.42 tor 
Ra* = 10 ^ and then increases to 0.43 for Ra* = 10 . This variation in the 
energy £atio is nearly negligible. Calculation (not presented here) a so 
shows Nuh/(4 Ar Nh) is equal to 0.56 for Ra* = 0 (heat conduction only). 

Thus, the energy ratio Nub/(4 Ar Nh) is significantly reduced by the buoyancy 
effect at lower modified Rayleigh numbers. The above analysis is obtained 
under the assumption of a flat interface. If the change of Ra is due to 
the change in gravity, the above qualitative conclusion may apply only when 
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the configurations of the interface are approximately the same, 
Bond number (Bo) is much greater than one. 


i .e . , the tank 


„ on In * he p F es ? nt calculations, the Prandtl number is varied from 1.25 to 
tt solutions (not presented in the paper) show that Pr has only a 
slight effect on the flow pattern and temperature distributions. A higher 

t enhances * be bu °yancy effect, resulting in greater radial veloc- 
lty at the interface. Thus, the maximum temperature, located at the tank 
sidewall closer to the interface, decreases with increasing Prandtl number. 
Following the above observation, the energy ratio, also seen from figure 4, 

tank 6 wfth f I Om - 0 1 ' 4 n f °A dV 1,2 r to 0 41 for Pr = 2 -65 in a uni formly-heated 
tranclWc ^ ~ J‘?i, and Ra = 106 ‘ Therefor e, even though an increase in Pr 
transfers more of the incoming energy into the near-interface region, the 
degree of superheat in the liquid is reduced. 


Effects of Nh and Ar 

inrrp«i!f!h ase }? £*“♦*?! 1 heat_flux Parameter, Nh, can be obtained either by 
increasing the wall heat flux, q w , or by decreasing the temperature difference, 

Us - *bl; Higher q w yields warmer fluid in the natural convection thermal 
? yCr T°T ing £° Ward . he interface - A lower value of temperature dif- 
!ho oK«’ T ?r~ ! b) : r ® duce f the heat transferred to the cold bottom. Both of 
the above effects tend to have more of the incoming energy transferred to the 

s?^if?£!^H Ce r T° n - Therefore ’ the energy ratio, Nu b /(4 Ar Nh), decreases 
gnificantl y with increasing Nh. Figure 5 shows the variation of the energy 
ratio as a function of the wall heat-flux parameter, Nh, for Ra* = 10 6 Also 

f?lro“ Pa p 1S ° n ’ - he pu £ e . conduct ion solution (i.e., Ra* = 0) is included in the 
gu . or a given Ra , if Nh is small enough, the buoyancy effect on the 
energy ratio is negligible. For Pr = 1.25, Ar = 1, and Ra* Jufi, the buoy- 
ancy has no effect on the energy ratio for Nh = l as seen from figure 5. This 
critical Nh decreases with increasing Ra* . 

the ® also yields smaller temperature difference between 

•*2? n d J he fluid ’ 1 - e -- the maximum temperature is decreasing 
U 8 Nh 'i ThlS results in low velocity buoyancy- induced boundary 
ayer flow and a weaker counter-clockwise vortex. Since the momentum of the 
urning flow at the centerline region is not strong enough to overcome the 
fraai ad y e f® e P ressu re gradient, a small separation bubble is formed in the 

*7 7 f j f f th ^ C T e l U ™- T1 ? e returnin g fIow moves radially toward the tank 
wall and feeds the boundary layer, resulting in the counter-clockwise vortex 

u PP e f Portion of the tank. Figures 6(a) and (b) shov, that 
the lower half of the tank is nearly a dead zone and is dominated by heat con- 
duction for a value of Nh = l. * 

The effect of the tank aspect ratio (Ar) on the flow pattern of a 
umlormly-heated tank is shown in figures 7(a) to (d) for Pr = 1.25, Nh = 4 
and Ra = 10 6 . The flow pattern for Ar = 1 is shown in figure 2(b). The 
counter clockwise vortex becomes more symmetric, with the vortex center moving 
from the upper right corner toward the midheight of the liquid region, as Ar 
is increasing. Also, an opposite weak vortex is formed at the centerline 
region and moves quickly toward the tank bottom as Ar increases. 
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Figures 8(a) to (d) show the effect of Ar on the temperature field. The 
maximum temperature located at the tank sidewall increases with Ar. Gener- 
ally, the liquid near the interface region is superheated over the entire tank 
cross-section. However, if Ar is too small, the convection boundary layer 
will not be well developed because the distance between the tank bottom and 
the interface is too short. Thus, only a small region confined in the upper 
right corner can be superheated. This phenomena is shown in figure 8(a) for 
Ar = 0.25. Figure 9 shows that a lower Ar results in greater value of Nu^/ 
(4 Ar Nh). It is noted that if Ar is too small (e.g., Ar = 0.25 in fig. 9), 
heat is coming in through the interface instead of going out. Thus, the value 
of the energy ratio, Nu^/ ( 4 Ar Nh) , can be greater than one. 


Effect of Rq 

In this study, calculations are also made to investigate the effect of the 
heat-flux distribution parameter, Rq = y^w.l’ on the flow pattern and tem- 

perature fields. Figure 10 shows a significant decrease in the energy ratio, 
Nub/ (4 Ar Nh), as Rq increases. It is because more heating on the upper por- 
tion of the tank wall (R q > 1) enhances the buoyancy effect and contributes 
more heat in the near interface region. If the lower portion of the tank wall 
has more heat flux coming in (Rq < 1), the fluid near the lower portion of the 
tank wall moves upward in the boundary layer and encounters a relatively low 
temperature region near the upper tank wall where the wall heat-flux is lower. 
This results in a mixing between the bulk liquid and the natural convection 
boundary layer. Thus, the heat transferred to the interface is reduced and 
Nub/(4 Ar Nh) is increased. From this observation, one can expect that if the 
given average heat-flux can be arranged such that Rq < 1, the tank pressure 
rise rate will be lower. This finding is also in agreement with the experimen- 
tal work by Evans et al. (ref. 7). 


Figures 11(a) and (b) show the streamline plots for Rq = 1/3 and 3, 
respectively. The streamline plot for Rq = 1 is shown in figure 7(c). As Rq 
increases from 1/3 to 3, the center of counter-clockwise vortex moves toward 
thq interface resulting in a higher radial interface velocity. The weak sec- 
ondary vortex tends to detatch from the tank bottom. The isotherm plots for 
Rq = 1/3 and 3 are shown in figures 12(a) and (b) , respectively. The plot for 
Rq = 1 is shown in figure 8(c). For Rq = 1/3, there are two peak temperatures 
at the tank wall with the lower peak being the maximum temperature. As Rq 
increases, the value of lower peak decreases and the value of the higher peak 
increases becoming the maximum temperature. If Rq is large enough, the lower 
peak disappears and the maximum temperature moves slightly upward. On the 
other hand, if R q is decreasing from 1/3, the value of upper peak is decreas- 
ing and the lower peak keeps increasing. The upper peak then disappears when 
Rq is small enough. 

At the interface, the numerical solutions show the local value of u* is 
generally two orders of magnitude smaller than v* (except in the small regions 

near the tank sidewall and centerline), especially for large Ar and Nh and 
low Ra*. Thus, it is justified that the neglect of u e in the calculations 
is acceptable. For the problems considered in the study, a large portion of 
liquid in the tank is superheated. Therefore, the above analyses are meaning- 
ful only if there is no boiling in the liquid, i.e., the maximum temperature, 
T max> is less than the incipient boiling temperature. A practical application 
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?h f J h i S .fK d V S t0 pr ® dict t 5 e required subcooiing for a given heat flux so 
that bulk boiling can be avoided. 


CONCLUSIONS 

r.art;In S ?*?? r P reaen * s . a numerical study of buoyancy- induced convection in a 
partially filled cylindrical tank subjected to wall heat-flux. The liquid- 
vapor interface is assumed flat, shear free, and at a constant saturation tem- 
perature. The tank bottom is kept at a constant temperature lower than the 
interface temperature Several parameters governing the fluid and thermal 
behavior are identified and their effects on the velocity and temperature 
fields are investigated The results of the investigation based on the param- 
eter ranges covered in the present study are summarized as follows: 

f 1 * F ? r t uniforra heat flux - the flow pattern usually contains a main 
^: er - Cl ? CkWi r. V0 ^ eX ca rT ng superheated liquid. Also, a weak clockwise 
vortex is formed in the tank bottom region near the centerline. In the tank 

bottom region, heat conduction is generally dominant because the fluid velocity 
is very I ow , 

• tv. 2 ' The ener sy ratio : Nu b /( 4 Ar Nh), is a fairly weak function of Ra* 

in the range considered in the study. The energy ratio decreases as Pr 
increases . 


3. Maximum liquid temperature occurs at the upper portion of the tank 
sidewall. For a given Nh and Ar, the maximum liquid temperature decreases 
with an increase in Ra and Pr. 

4. An increase in Nh and Ar results in an increase in the value of 

rwlnlfn liq +i! d J em P e ^ ature aad a decrease in the energy ratio, Nub/(4 Ar Nh). 
Generally, the liquid near the interface region is superheated over the entire 

fnL Cr ° S ^ SeCtl0n - , H ° wever - if Ar is too small, only a small region con- 
in e u PP er right corner can be superheated and the energy ratio may be 
greater than one. * 

5 . If Rq is large enough, the maximum liquid temperature is located at 
the upper porfion of the tank sidewall. As R q decreases, the value of the 
maximum temperature decreases while a second temperature peak occurs at the 
lower Portion of the tank sidewall. Further decrease in R a , the upper temper- 
ature peak disappears and the lower temperature peak increases and becomes the 
maximum temperature of the liquid. As a result, the value of Nuk/( 4 Ar Nh) 

increases with decreasing Rq due to the mixing between the convection bound- 
ary layer and bulk liquid. 
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Figure 1 . — Physical system and coordinates. 
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Figure 2. — Streamline plot for various Ra*. (Ar = 1.0, Pr = 


(a) Ra* - 10 5 . (b) Ra* - 10 6 . 



Figure 3.— Isotherm plot for various Ra*. (Ar - 1 .0, Pr - 1.1 
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(c) Ra* « 1 0 9 . 
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Figure 4 —Effect of Rayleigh Number and Prandti Number 
on the fraction of the incoming energy transferred through 
the cold surface. 



Figure 5— Effect of wall heat-flux parameter on the fraction of 
the incoming energy transferred through the cold surface. 
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Figure 6. — Streamline and isotherm plot for 


(b) Isotherm. 



= 1 , Ar = 1. Ra* = 10 6 , Pr = 1.25, and Nh = 1. 
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Figure 9.-— Effect of aspect ratio on the fraction of the in- 
coming energy transferred through the cold surface. 
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Figure 10. — Effect of wall heat-flux distribution parameter 
on the tract incoming energy transferred through cold 
surface. 
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(a) Rq = 1/3. 



Radial distance, r* 


(b) Rq = 3 



Figure 1 1 .—Streamline plot for Rq = 1/3, and 3. (Ar = 0.75, Ra* = 10 6 , Pr = 1 .25, and Nh = 4) 


(a) Rq = 1/3. (b) Rq = 3 




Figure 12.— Isotherm plot for Rq = 1/3, and 3. (Ar = 0.75, Ra* =10 6 , Pr = 1.25, and Nh = 4) 
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